Mechanisms of two-color laser-induced field-free molecular orientation 
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Two mechanisms of two-color (lu + 2lj) laser-induced field-free molecular orientation, based on the 
hyperpolarizability and ionization depletion, are explored and compared. The CO molecule is used 
as a computational example. While the hyperpolarizability mechanism generates small amounts of 
orientation at intensities below the ionization threshold, ionization depletion quickly becomes the 
dominant mechanism as soon as ionizing intensities are reached. Only the ionization mechanism 
leads to substantial orientation (e.g. on the order of {cos9) > 0.1). For intensities typical of 
laser-induced molecular alignment and orientation experiments, the two mechanism lead to robust, 
characteristic timings of the field-free orientation wave-packet revivals relative to the the alignment 
revivals and the revival time. The revival timings can be used to detect the active orientation 
mechanism experimentally. 



Laser-induced field-free alignment of small gas phase 
molecules, where the molecular axis is aligned along par- 
ticular direction (see FigH]), is now routine [l[ and is a 
quickly bccominga central tool to in attosecond f3| and 
photoionization [3| experiments. Much less studied is 
the laser-induced field-free orientation of polar molecules, 
where both the molecular axis and the asymmetry point 
along a particular direction, which has recently been 
achieved 3 using two-color ultrafast pulses built from 
a fundamental and its second harmonic. The under- 
lying physical effect thought to be responsible for the 
orientation is the hyperpolarizability interaction 0, 
However, it is possible that a second mechanism - ion- 
ization depletion - is operative, where a two-color ul- 
trafast pulse selectively ionizes molecules at particular 
angles with respect to the polarization direction of the 
laser field. Hyperpolarizability generates orientation by 
causing an asymmetrical force that pushes the molecules 
toward orientation, while ionization depletion generates 
orientation by directly carving out an asymmetrical an- 
gular distribution upon ionization. In this work we im- 
plement a simple model of the ionization depletion mech- 
anism, and contrast and compare the two mechanisms. 
A rigid linear rotor is used as an example system. Exten- 
sion to symmetric and asymmetric tops is conceptually 
straightforward. Rotor parameters corresponding to the 
CO molecule are used in numerical examples. 

The two-color laser pulse is written as 



E{t) = EQf{t)[cos{ujt) + cos(2a;i)], 



(1) 



where /(t) is the pulse envelope (0 < f{t) < 1), and Eq 
(a positive real number) controls the peak electric field 
strength. The relative phase of the two fields was set 
to zero, and their relative amplitudes were set to one, 
in order to maximize the field asymmetry that leads to 
orientation. The nuclear rotation effectively does not re- 
spond on the timescale of the carrier oscillations, and it 
is then appropriate to use the cycle-averaged rotational 



Hamiltonian of the system: 

H{e,t) = BJ{J+l) + Vp{0,t) + VH{e,t) + Vi{6,t), (2) 

where B is the rotational constant, Vp{9,t) is the po- 
larizability term that generates molecular alignment [Q], 
VH{d,t) is the hyperpolarizability term, and Vi{9,t) ac- 
counts for ionization (see below). All equations use 
Hartree atomic units (me = e = h — 1). For E{t) in 
Eq. [Tl the two middle terms in Eq. [5] are given by 5] 
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where Aa = a|| —a± is the polarizability anisotropy, and 
the Pijk are elements of the hyperpolarizability tensor. 
The molecular constants chosen to model CO are given 
in Table m 

Ionization is introduced using a complex absorbing po- 
tential 



Vii9,t)^-{i/2)Ti{9,t) 



(5) 
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FIG. 1: Cartoon illustrating the difference between alignment 
and orientation of a schematic polar diatomic molecule. 
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TABLE I: Molecular constants (a.u.) used to model CO. 



Parameter 


Value [Ref.] 


Parameter 


Value [Ref.] 


B 


8.7997x10"'' [3 


Aa 


3.6 [8] 




28.91 [81 


Pxxz 


7.69 [8j 


Ip 


0.516 [7] 


Co 


0.2214 xlO"^ 


Cl 


-0.2141 xlO"^ 


C2 


0.0822 xlO~^ 
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where TjiO^t) is the cycle-averaged ionization rate. The 
remaining time-dependence of Ti[6,t) arises from the 
pulse envelope f{t). The complex potential causes non- 
unitary quantum evolution that removes amplitude as a 
function of angle, modeling the effects of population loss 
due to ionization. For simplicity, a separable form of 
Ti{6, t) is assumed 
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Vi(e,t)^~{i/2)K{t)T,,i(e), 



(6) 



where rref(^) is the cycle- averaged angular ionization 
rate calculated at a characteristic field strength i?rcf , and 



K{t) = exp --(2/p)3/2 [\Eofit)\-' - \EU 



(7) 



is the tunneling exponent 0, [lol | that provides the dom- 
inant scaling of strong field ionization. Ip in Eq. [S] 
is the ionization potential of the molecule. Scaling in 
Eq. [6] qualitatively captures the ionization yield in small 
molecules over several orders of magnitude in laser inten- 
sity [Hi]. 

Calculation of Tjciid) for ionization of X^I]+ CO to 
X2S+ C0+ is carried out using the time-dependent 
mixed orbital/grid method of Ref. in the single- 

channel approximation. The wavefunctions of the neutral 
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FIG. 2: Cycle-averaged ionization rate Frcf (6') Shown arc the 
numerically calculated rates, and the fit from Eq. 
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FIG. 3: Left axis: the maximum field-free orientation 
(j(cos6')j) after application of the two-color pump, as a func- 
tion of pump intensity. Each of the two colors in the pump 
pulse has intensity Jq. Thin line: the hyperpolarizability 
mechanism. Thick line: ionization depletion. Right axis: to- 
tal ionization probability (dashed line). The vertical dotted 
lines gives the intensities used in Figs. |3] and [S] 

and the ion are calculated using GAMESS-US fis!] with 
the cc-pVTZ basis set [1^] at the complete active space 
level using 10/9 (neutral/cation) active electrons in 8 or- 
bitals. Uniform Cartesian grid extended to ±13 Bohr, 
with spacing of 0.2 Bohr. The time step was 0.002666 
a.u. The simulation was run for a full cycle of the two- 
color field Eq. dH), with Eo = = 0.0535 (Jo = lO^* 
W/cm2), w = 0.057 (A =800/400 nm) followed by 2 fs at 
zero field. The cycle-averaged ionization rate was calcu- 
lated by integrating the flux absorbed [15| at the edges 
of the grid, and dividing by the cycle duration (2tt/u}). 
Fig. m plots the calculated T^ci{6). Also shown is the fit 
given by a truncated Fourier series (see Table U]): 

rrcf(6') = Co +C1COS0 + C2COS26' (8) 

All subsequent calculations use the fit of Eq. |5] in Eq. [51 
Although the Hamiltonian ([2]) includes both the hy- 
perpolarizability and ionization terms, numerical results 
below consider the two mechanisms individually, to elu- 
cidate the characteristic features in the induced orienta- 
tion. Simulations with both mechanisms active simul- 
taneously did not reveal any qualitatively new features, 
and will not be discussed. 

The initial populations of rotational states \J,M) are 
given by the Boltzmann distribution at temperature T — 
50K: 

exp(-BJ(J + l)/fcT) 



PiJ,M) 



(9) 



Ej'C^J' + l)exp(-BJ'(J' -t- l)/kT) 

The time evolution of each rotational state within the 
ensemble is expanded in a spherical harmonics basis 



|*(<))=^a,;W|J,M). 



(10) 
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FIG. 4: The hyperpolarizability orientation mechanism. 
Alignment (thin line, left axis) and orientation (thick line, 
right axis) wave-packet revivals for a selection of intensities. 
The vertical dashed line denotes the revival time trev 



The single sum over J is appropriate since the Hamilto- 
nian in Eq. ([2]) conserves M. The Schrodinger equation 
for the coefficients a,j{t) is 



= BJ{J+l)aj{t) 

ot 



(11) 



J2{J, M\Vp + Vh + Vi\J\M)aj,{t). 



The envelope function f{t) is defined as 



fit) 



0, 

sin (7rt/2To. 
0, 



t < 

0<t < 2Ton 



(12) 



corresponding to a sin^ pulse for the intensity I ~ E'^. 
The parameter Ton {tou = 30 fs) is the full width at half- 
intensity. While the pulse is on (< < 2ro„), Eq. (fTTj) 
is solved using the Crank-Nicholson method [16]. After 
the end of the pulse {t > 2to„), the analytical field- free 
propagation solution is used. Each initial \J,M) state is 
propagated independently. The observables are averaged 
over the thermal distribution P{J,M). 

The maximum field-free orientation following the pulse 
is shown in Figure [3l as a function of the single-color in- 
tensity /o — Eq. The orientation is characterized using 
the observable (cos 6) . The total ionization probability is 
also shown (dashed line, right axis) on the same figure. 
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FIG. 5: The ionization depletion mechanism. See Fig. U 



The peak orientation scales very differently with the in- 
tensity for the two mechanisms. The hyperpolarizability 

3/2 

mechanism shows the expected /q scaling. The deple- 
tion mechanism shows the same scaling as the ionization 
probability, coming from the tunneling exponent (Eq. [7]) . 

At low intensities, where ionization is not possible hy- 
perpolarizability provides the dominant mechanism of 
orientation generation. When the intensity increases and 
substantial ionization occurs, ionization depletion domi- 
nates the orientation. For the present model, the cross- 
over occurs at 1.05 x 10"'^* W/cm^. This cross-over inten- 
sity is sensitive to the exact molecular parameters used, 
especially to the ionization rate model. Although quali- 
tatively reasonable, Eqs. ([6])-([8]) are not expected to be 
quantitatively accurate, and consequently gives only a 
rough estimate of the experimental cross-over intensity. 
Due to the rapid increase in the peak orientation once 
ionization is active, it is likely that ionization depletion 
will be the active mechanism when the largest degrees of 
orientation are reached, regardless of the intensity where 
the mechanisms switch. Further, since the peak orien- 
tation for the ionization mechanism parallels closely the 
ionization probability, it follows that substantial orienta- 
tion requires substantial ionization. For example, reach- 
ing an orientation of {cos0) ^0.1 requires ionizing more 
than 10% of the sample in the CO case. 

During short-pulse alignment and orientation, the 
molecular ensemble undergoes wave-packet revival dy- 
namics, where short periods of sharp angular localization 
leading to strong alignment and orientation (the wave- 
packet revivals) are separated by longer periods charac- 
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terized by dispersed wavefunctions. Figures |4] and [5] show 
the ahgnment and orientation revivals at the first full re- 
vival trev = tt/B ss 8.64 ps. The alignment is character- 
ized by the commonly-used alignment parameter (cos^ 0) . 
The orientation is tracked using (cos 6*)^, the square of 
the directional orientation parameter {cos 9). Note that 
{cos 9)^ does not contain the orientation direction. This 
is the appropriate observable to characterize some exper- 
imental orientation probes, such as high-harmonic gen- 
eration For the present study, this observable con- 
veniently draws out the peak orientation independent of 
direction. 

At all intensities, the hyperpolarizability mechanism 
leads to peak orientation at early times relative to trev 
(Fig.|4]). Since Uev is also typically near the center of the 
alignment revival, the shift of the orientation wavepacket 
relative to trev can also be seen as a shift toward, and be- 
yond, the minimum of the alignment revival. Ionization 
depletion, on the other hand, leads to the peak orien- 
tation shifted to later times relative to trev and to the 
alignment revival (Fig. [S]). Thus, the relative timing be- 
tween the peak of orientation and the alignment revival 
can be used to determine the active mechanism of orien- 
tation experimentally. 

It is instructive to examine the origin of the relative 
timings within the two orientation mechanisms analyt- 
ically, using the impulsive approximation. If the pulse 
duration is much shorter than typical rotation timescales 
{Ton ^ ^/B), the Hamiltonian of Eq. [2] becomes 

H = BJ{J + 1) + [F/.(0) + Vlf{9) + Vlie)] 6{t). (13) 

where VJ{6) ~ J Vi{9,t)dt. Without sacrificing any of 
the essential physics of the model, the potentials V- can 
be approximated by the leading non-constant terms in 
the cosO expansion 



Vl,i9) 

VR0) 



-Xp cos^ ( 
-Xh cos 6 
-iXi cos 6 



(14a) 
(14b) 
(14c) 



where the numerical parameters (at the single-color in- 
tensity Iq = 6.0 X 10^3 w/cm^) are Xp = 3.823, = 
0.2561, and A/ = 1.095 x lO'^. The t = +0 wavefunction 
then becomes [HI 

*) = exp [iXp cos^ 9 + {iXh - A/) cos 6*] |*o)- (15) 

Impulsive orientation revivals for Iq — 6.0 x lO^'^ W/cm^ 
and the initial ensemble temperature of 50K are shown 
in Fig. [6l The impulsive approximation reproduces the 
temporal structure of the revivals seen in the full sim- 
ulations (cf. the bottom panel of Fig. [6] to the bottom 
panels of Figs. H] and [5|). Thus, the temporal signature of 
the orientation is unrelated to the finite pulse duration 
effects. Furthermore, the orientation revival timings are 
not affected by the fine details of the angular dependence 
of the field- molecule interaction, and are robust. 
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FIG. 6: Orientation revivals within the impulsive approxima- 
tion (Jo ~ 6. Ox 10^"^ W/cm^). The ionization depletion results 
are scaled by a factor of 1000. Top panel: revival structure 
with the alignment interaction turned off (Ap = 0). Bottom 
panel: with the alignment interaction included. The vertical 
dashed line denotes the revival time trev 



Consider now the orientation generated with align- 
ment interaction turned off (Fig. |6l top). Both mech- 
anisms generate revivals that are completely symmetric 
around trev ■ The hyperpolarizability revivals have a zero 
at trev and peak a short time after. The ionization deple- 
tion generates an orientation peak exactly at trev The 
alignment-free behavior of the orientation is easy to un- 
derstand. The rotational wavepacket is periodic, with pe- 
riod of trev Thus bchavior at trev is a 'mirror' of the ini- 
tial motion. Within the hyperpolarizability mechanism, 
the molecules get an instantaneous kick from the S(t) 
pulse, but the spatial orientation distribution remains 
unchanged. It takes a short time for the imparted angu- 
lar momentum to cause the molecules to rotate toward a 
point of maximum orientation. This is the exact behavior 
seen in Fig. [5] Within the ionization depletion mecha- 
nism, the S{t) pulse instantaneously generates non-zero 
orientation by directly modulating the angular distribu- 
tion. The instantaneously-generated orientation should 
then disperse as the molecules start to rotate. Again, 
this is clearly seen in Fig. [6] 

Including now the effects of alignment to the laser- 
induced orientation results in the revival structures seen 
in the bottom panel of Fig. [S] The hyperpolarizability 
revival becomes asymmetric, but does not change posi- 
tion, while the ionization depletion revival keeps its basic 
shape, but shifts to later times. The coupling of the align- 
ment to the ionization depletion revivals is qualitatively 
clear. The S{t) pulse first creates instantaneous orienta- 
tion. Then, as the cos^ 9 interaction generates a torque 
that squeezes the angular distribution toward molecular 
alignment, the instantaneously-generated orientation is 
also squeezed toward greater orientation as the alignment 
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maximizes (Figs. |4l [5]). The coupling between the align- 
ment dynamics and the hyperpolarizability-generated 
orientation revival is not as obvious, but it is still ac- 
counted for by the coupling between the aligning and 
orienting forces. 

Two mechanisms of laser-induced molecular orienta- 
tion, due to hyperpolarizability and ionization depletion, 
were studied using the CO molecule as a model sys- 
tem. At low intensities, the hyperpolarizability mech- 
anism dominates the generation of molecular orienta- 
tion. Once the ionizing intensities are reached, the ion- 
ization depletion mechanism dominates. The two orien- 
tation mechanisms have clear parallels to the well-known 
sources of vibrational wave packets: impulsive stimulated 



Raman scattering (ISRS) [l9| . and R-dependent ioniza- 
tion ("lochfrass" [23). Indeed, the hyperpolarizability 
mechanism acts by imparting angular momentum on the 
system, similar to ISRS. The ionization depletion on 
the other hand, directly "burns-in" the rotational wave 
packet in the coordinate representation, analogously to 
lochfrass. The orientation revivals contain an unambigu- 
ous signature of the active orientation mechanism, the 
hyperpolarizability mechanism, the peak of the orienta- 
tion revival precedes the revival time ire-u, while for the 
ionization depletion mechanism the peak of the orienta- 
tion revival lags behind treu- These signatures arise due 
to interplay of the aligning and orienting interactions. 
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